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ELASTIC  SOLUTIONS  IN  A  SEMI-INFINITE  SOLID 
WITH  AN  ELLIPSOIDAL  INCLUSION 

1.  INTRODUCTION 

The  elastic  fields  due  to  inclusions  in  infinite  media  have  been  extensively  investigated 
using  Eshelby’s  method  (  1957,  1959,  1961  ).  Cases  of  more  practical  interest  are  the 
elastic  fields  due  to  inclusion  in  semi-infinite  media  since  many  sites  of  initial  strains, 
resulting  from  high  stresses  due  to  contact,  heating,  or  metallurgical  transformations,  are 
confined  to  a  near  surface  zone.  The  problem  of  an  ellipsoidal  inclusion  which  has 
undergone  a  simple  shear  is  of  interest  in  connection  with  twinning  and  martensitic  or  other 
diffusionless  transformations.  The  elastic  solution  for  an  inclusion  near  the  free  surface  has 
been  solved  for  a  spherical  inclusion  with  pure  dilatational  eigenstrain  (  stress  free 
transformation  strain  )  (  Mindlin  and  Cheng,  1950B  ),  an  ellipsoidal  inclusion  with  pure 
dilatational  eigenstrains  (  Seo  and  Mura,  1979  )  and  a  cuboidal  inclusion  with  uniform 
eigenstrains  (  Chiu,  1977  ).  In  both  Mindlin  and  Cheng's,  and  Seo  and  Mura's  analyses, 
the  solutions  are  obtained  by  integrating  the  Green's  function  of  a  point  force  in  the  interior 
of  a  semi-infinite  solid  (  Mindlin,  1936,  1953  ).  In  Chiu’s  analysis,  the  free  surface 
condition  is  satisfied  by  superimposing  the  solution  of  a  half-space  under  normal  surface 
traction  on  the  full  space  solution  due  to  a  cuboidal  inclusion  and  its  image  with  the  uniform 
eigenstrains. 

.The  solution  of  the  linear  equations  of  equilibrium  of  an  elastic  body  with  a  force 
acting  at  a  point  within  an  isotropic  body  bounded  by  a  plane  has  been  solved  by  starting 
with  Kelvin's  solution  for  a  force  in  an  infinite  body  and  guessing  the  nuclei  of  strain  to 
add  outside  of  the  semi-infinite  body  so  as  to  annul  the  tractions  on  the  plane  boundary  ( 
Mindlin,  1936  ).  The  same  results  have  been  obtained  directi'  V  means  of  an  application 
of  potential  theory  (  Mindlin,  1953  ).  The  stress  functions  of  dii  t  strain  of  nuclei  in  the 
semi- infinite  elastic  solid  (  Mindlin  and  Cheng,  1950A  )  are  derived,  by  the  processes  of 
superposition,  differentiation  and  integration,  from  the  solution  for  the  single  force  in  the 
interior  of  the  semi-infinite  solid. 

Maauacript  approved  November  13,  1989. 
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The  method  of  Hankel  transformations,  elaborated  for  cylindricaily  symmetrical 
problems  of  the  theory  of  elasticity  in  Sneddon's  book  ( 1951  ),  has  been  used  to  solve  the 
stress  field  of  a  circular  edge  dislocation  loop  with  Burger's  vector  normal  to  the  plane  of 
the  loop  ( prismatic  loop  )  in  an  unbounded  solid  (  Kroupa,  1960  )  and  in  the  half  space  ( 
Bastecka,  1964 ). 

In  the  present  study,  Eshelby's  method  for  the  ellipsoidal  inclusion,  the  Hankel 

transformation  method  for  the  axisymmetric  problems  and  the  Mindlin’s  stress  functions 

for  strain  nuclei  of  double  force  with  moment  are  used  for  the  analysis  of  the  elastic 

solution  of  an  ellipsoidal  inclusion  in  the  half  space  when  a  uniform  eigenstrain  with 
T  T  T  T  T  T 

components  e33 ,  etl=  e22,  e31 ,  e23  and  e12=  0  are  given  initially  inside  the  inclusion. 
Existing  solutions  are  shown  to  be  special  cases  of  the  present  one. 


2.  ELASTIC  SOLUTIONS 

T 

The  present  problem  is  to  express  the  elastic  field  when  the  eigenstrain  e-  in  an 

ellipsoidal  subdomain  (with  semi-axes  at ,  a^,  a^,  and  center  at  x,  =  x2  =  0  and  x3  =  c  ) 

T  T  T 

of  the  half  space  x3  >  0  (  Fig.l  )  ts  made  up  of  components  e11  =  e22,e33, 
T  T  T 

e3 1 ,  e23  and  e  l2  =  0.  In  order  that  the  plane  x3  =  0  be  a  surface  free  of  external  tractions, 
the  stress  components  on  this  plane  must  satisfy  the  following  boundary  conditions 

(  <J13)x,  =  0=  (  **23  )*3=  0  =  (<y33)x3  =  0=®»  (1) 

and  both  the  equations  of  equilibrium 

%.r  °- 

and  the  compatibility  equations 
2  1 

V  o..  + - aH  =  0,  (3) 

J  1+v  ,IJ 


2 


where  the  numerical  suffixes,  i,  j  =  1,  2,  3,  following  a  comma  denote  differentiation  with 
respect  to  the  Cartesian  coordinates  xlf  x2,  x3,  respectively;  a  repeated  suffix  is  summed 
over  values  1,  2,  3,  and  v  is  the  Poisson's  ratio. 

Similar  to  the  approaches  ofBastecka  (  1964  )  and  Chiu  (1978),  the  stress  in  the  half 
space  x3  >  0  outside  the  axisymmetric  ellipsoidal  inclusion  centered  at  the  point  (0,0,c)  can 


be  expressed  as 


i  n 

a.  ■  =  o. .  +  <y.  •  +  a ■  ■ , 
»j  »j  ij  »j’ 


(4) 


which  satisfies  the  required  boundary  conditions,  Eq.(l),  the  equilibrium  condition, 
Eq.(2),  the  compatibility  equations,  Eq.(3),  and  also  converges  to  zero  for  x,  and  x2 


l  . 


approaching  ±oo  and  x3  approaching  oo  .  In  Eq.(4),  the  term  |  is  the  stress  caused  by  the 
inclusion  centered  at  (0,0, c)  in  an  isotropic  infinite  body  and  can  be  obtained  by  using 
Eshelby's  method  (1961)  with  displacements 


„■= 


1  T  I  1  T  I  V 

8tt(l-v) ejkV  ,jk  2k  C‘  ,k  4k(1-v) 


Tl  I 


(5) 


i  i 

where  V  and  9  ,  respectively,  are  the  biharmonic  and  harmonic  potentials  of  attracting 

T  T 

matter  of  unit  density  filling  the  volume  Qj  and  e  =  emm  with  repeat  index  m  sum  over  1 ,  2 

n 

and  3.  is  the  stress  caused  by  the  image  inclusion  Gt^  centered  at  the  point  (0,0,-c)  in  an 
isotropic  infinite  body,  with  eigenstrain 


(e^)n=-  (e-)1.  (6) 

The  solution  for  the  stresses  a  °  is  obtained  by  translating  the  origin  of  coordinates  in 
Eq.(4)  from  points  (0,0, c)  to  (0,0,-  c).  The  additional  stress  in  Eq.(4)  is  the  fictitious 
stress  necessary  to  make  the  surface  of  the  half-space  free  of  stresses  and  it  satisfies  the 
boundary  conditions 


(  °  33)*!  =°“  ( 


a33+  <733)x, 


3 


(7) 


(  <7  13  =0 

( ®  23^xi  =  0 

The  elastic  solution  of  an  inclusion  in  a  half-space  present  in  this  report  are  solved  first 

T  T  T 

by  considering  the  inclusion  with  pure  principal  eigenstrains  ej  j  =  e22  >  e3  3  then  the 

T  T  T 

inclusion  with  pure  shear  eigenstrains  e31 ,  e23  and  e12  =  0  is  solved.  The  final  solution  is 
the  linear  superposition  of  these  solutions. 

T  T  T 

(A)  Inclusion  with  with  pure  principal  eigenstrams  e  1 1  =  e22 ,  e33 . 

The  stress  field  of  the  axisymmetric  penny  shape  inclusion  with  principal  eigenstrains 
T 

e33  in  an  infinite  medium  obtained  by  Eshelby  method  is  compared  with  the  stress  field  of 
a  prismatic  loop  with  radius  a  and  Burger's  vector  b  in  an  infinite  medium  as  obtained  by 
Kroupa  (  Yu  and  Sunday,  1988  ).  A  relationship  is  found  between  the  potential  function  <j> 
of  the  inclusion  and  the  integral  function  I0\  which  involves  the  product  of  Bessel 
functions  Jm ,  for  the  solution  of  the  prismatic  loop.  That  is 

Iq  = - e33<l,»  (8) 

2tiab 

where  in  cylindrical  coordinate  (  r,  9,  z  ) 

.  (9) 

and  Jm  is  the  Bessel  function  of  the  mth  order.  In  Eqs.  (8)  and  (9),  both  the  harmonic 
potential  $  and  the  integral  function  Iq1  are  taken  the  origin  as  the  center  of  the  inclusion 
and  the  center  of  the  dislocation  loop,  respectively. 

Then  the  fictitious  stress  field  is  solved  first  for  the  two  dimensional  problem  by  using 
the  Hankel  transformation  method  and  then  it  is  transformed  into  the  three  dimensional 
case  by  use  of  the  relationship  between  and  Iq1.  The  displacements  solved  are 


.  1  a  X 

—  ~  \  ^i3)xj=0» 

/  1  n  \ 

—  a23+ 


4 


,  T  T 

v  fin  j  i  i  n  n 

U1  =  [  ¥,133-  2vfl+  2zV.1333  +  (3-4v)V.!33 


8tc(1-v) 


-  2z2<J>^33_  4(2-v)z<(>  ^ 3  —  2(2— 3v)4>  ^  ] 


(l+v)en  t  n  n 

- [<j>  i+2z(j)  13  +  (3-4v)  <J>  |  ], 

47t(l-v) 


T  T 

( c*»3  Ci  i )  j  j  n  n 

U2=  _____  [  ¥  233-  2v<|»>2+  2zV|/i2333+  (3-4v)Vi233 


8k(1-v) 


-  -  4(2-v)z't».23  -  2(2-3v)<j)  “  ] 


(l+v)eu  !  n  n 

- [<t>  2+2z^  23+ (3-4v)  ^,2]’ 

4it(l-v) 


T  T 

i  e,,-  c, , ;  i  I  it  n 

u3=  — -  —  [  V,333~  2(2-V)<|»  3+  2zy  3333-  (3—4v)¥  333 


87t(l-v) 


-  2z2<1>  333- 4(l+v)z<|»  33  +  2(4-5v)*“] 


(l+v)cu  j  n  n, 

[<j>  3+2z(|>  33-(3-4v)<|>  3], 


(10) 


4jt(l-v) 

and  the  stresses  are 

_  T 

m  e33  ®ii)i  i  n  _  a  11  n 

<*n= - [V,1133+2vf22+2zV  11333+ (3“4v)¥  1133- 4vV.3333 


C22  = 


4n(l-v) 

"  2z2<j>“133-  4(2-V)z(t>“  13  +  4vz(()  333-  2(2-3 v)  4»“  ,  +  14v  <t>^3  ] 
T 

(l+v)neu  in  n  n 

- [  <j>  11+  2z<|>  113+  (3-4v)  <{>  n-  4v  <t>  33], 

27c(1-v) 

T  T 

MX  e33  e  1 1 )  r  I 


4n(l-v) 


[  ¥  2233  +  2V(j»  ,  J  +  2zy  22333  +  (3^V)\|f  2233  “  4vV.3333 


-  2z2<j>  2233“  4{2-v)z(| >i223  +  4vz<1>,333-  2(2-3 v)  <t>  22  +  14v  $  33  ] 


(l+vjjie.i  in  n  n 

- [  $  22+  ^  223+  (3-4v)  $.22“  4v  ^,331’ 

2n(l-v) 
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<*33” 


T  T 

e33  eu)  i 

tV,3333 


4«(l-v) 


I  II  u 

-4f  33^.3333  +  4* 


,33 


n  2  n  n 

+  2zy  33333-  2z  <^.3333  ®z<^,333l 


<*12  = 


(l+v)|ie..  1  n  „  n 

[<j>  33“  <(>  33+  2z<|>  333  j, 

2ji(1-v) 

WeI3-eu)rJ 


4jc(1-v) 


tV.1233  “  2v<^.12+  ^zV,12333  +  (3-4v)Y,  1233 
—  2z  ^,1233“  4(2— v)  z^123  ~  2(2— 3v)  <t>  121 


(l+v)uei  1  in  n 

- [<)>  12  +  2 zij>  123  + (3-4v)(j>  12], 

2re(l-v) 


<*23  = 


^el3-eli)r  I 


4ji(1-v) 


tV,2333  ~  ^,23  +  2z\(f  23333  +  Y.2333 
-2z2<j>  2333-  8z4».233  “  2^.231 


(l+v)ncn  I  n  n 

— — - -  [  <t>.23+  2z<|>  233+  <j>  23  ], 

2tc(1-v) 


<*31  = 


^ela-eL)! 


1 


n 


4jc(1-v) 


[Y  ,1333  -  2$  13+  2z\y  13333+  Y ,1333 


2  n  n  _  n  „ 

-2z  ^,1333“  ^,133  ”2<|>  13] 


(l+v)pellrJ  ^  ^  i 

— —  —  [  <j>  13+  2z<|»  133+<j>  I3], 

2tc(1-v) 


and  the  dilatadonal  stress  is 

(l+v)fi(e33-ell)  1  n  _  n  _  ,  n  . 

o- - [  -  $  33+  7d>  33-  2\|f  3333  +  2z<|>  333  ] 


2ti(1-v) 
2(l+v)  Melt  " 


Jt(l-V) 


',33' 


(ID 


(12) 


where  y  °and  <^n  are  the  biharmonic  and  harmonic  potential  of  attracting  matter  of  unit 
density  filling  the  volume  Qj,  respectively. 
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The  ad  onal  justification  for  the  substitution  of  Eq.  (8)  is  that  Esheiby  (  1961  )  has 

shown  that  the  remote  field  of  a  finite  prismatic  loop  with  area  A  and  Burgers  vector  b  is 

the  same  as  the  remote  field  of  an  inclusion  of  arbitrary  shape  whose  volume  V  and  the 
T 

stress  free  strain  e33  ( parallel  to  b  )  satisfy 

Ve  J3  =  Ab,  (13) 

For  a  circular  edge  dislocation  loop  of  radius  a  and  the  x3-axis  (or  z-axis)  as  the  axis  of 
symmetry  in  an  unbounded  medium,  the  stress  field  is  found  by  Kroupa  (1960)  by  using 
Hankel  transformations.  Since  for  the  remote  field  the  inclusion  can  be  of  any  shape, 
therefore,  the  harmonic  potential  <j>  for  a  spherical  inclusion  of  radius  a  is  chosen.  Then 
Eqs.  (8)  and  (13)  give 

IqI=  fg;  (a->0),  (14) 

which  is  the  same  as  obtained  from  the  mathematics  formula  given  by  Eason,  Noble  and 
Sneddon  (  1955  ).  There,  the  substituting  of  Eq.(8)  is  not  limited  to  penny  shape  inclusion, 
it  can  apply  to  any  ellipsoidal  inclusion  because  we  can  let  the  radius  of  the  dislocation  loop 
approach  zero  and  integrate  the  solutions  over  the  volume  of  the  inclusion  which  gives  the 
results  in  term  of  potential  functions. 

.  T  T  T 

(B)  Inclusion  with  with  pure  shear  eigenstrains  e3  y ,  e23  and  e12  =  0. 

The  linear  elastic  solution  for  an  inclusion  with  pure  shear  eigenstrain  is  obtained 
indirectly  by  the  application  of  Kelvin's  solution  for  double  force  with  moment  in  an 
infinite  body  and  Mindlin's  solution  for  double  force  with  moment  in  an  semi-infinite  body 
.  A  comparison  between  Kelvin's  solution  for  double  force  with  moment  in  an  infinite 
body  and  Eshelby's  solution  for  an  inclusion  with  pure  shear  eigenstrain  is  conducted 


7 


first.  A  relation  between  these  two  solutions  is  established  first.  Then  the  solution  for  the 
inclusion  in  the  half  space  can  be  obtained  by  substituting  this  relationship  into  the  solution 
obtained  by  Mindlin  for  the  semi-infinite  solid. 


For  the  double  force  with  strength  A  in  the  x  {x  direction  with  moment  about  x2  axis 
and  the  double  force  with  strength  A  in  x  3  direction  with  moment  about  y  axis  in  infinite 
solid,  the  Galerkin  vector  is  ( Mindlin,  1936 ) 

F  =  iAx yR  +  kAx  j/R  ,  (15) 


2  2  2  2 

where  R  =  x  { +  x2  +  x3.  The  displacements  derived  from  Eq.  (15)  are 


A 

ut  =  — [2(l-v)<p  3-x3<p  u] , 

p 

u2  =  — [x3<p  12| , 

M- 

A 

u3  =  —  [  2(l-v)tp  13] , 

p 


where  <p  =  1/R. 


(16) 


For  penny  shape  inclusion  (disc;  a  t  =  a2,  a3  — >0  )  with  e3  { =  e  1 3  are  the  only  non- zero 
T 

components  of  e^.  Eq.  (1)  gives  the  displacement  as 


u  i =  ~ 

u2  =  - 


u3  =  - 


T 
e  13 

4tt(1-v) 

T 

e13 

4rc(l-v) 

T 

e13 

4n(l-v) 


[2(l-v)<j>  3-x3(|>  n] , 


I  X3^,12^  » 


[  2(l-v)<p1-x3<t»13] , 


(17) 


where  <)>  is  the  harmonic  potential  of  attracting  matter  of  unit  density  filling  the  inclusion 
volume  Q  centered  at  ( 0,0,0  )  and  is 
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$  =  JQ<pdT.  (18) 

Eq.  (18)  is  valid  for  all  shapes  of  ellipsoidal  inclusion  not  just  for  penny  shape  inclusion. 

By  comparing  Eqs.  (16)  "nd  (17),  it  is  found  that  the  solution  for  the  inclusion  with  pure 
T  T 

shear  eigenstrain  e31 =  el3  can  be  obtained  by  integrating  the  results  of  the  double  forces 
with  the  same  strength  A  in  both  the  x  t  and  x3  direction  with  moment  about  x2  axis  in  an 
infinite  solid  for  the  strain  over  the  volume  element  of  the  inclusion  provided  that 


A  =  - 


4n(l— v) 


For  the  double  force  at  point  (  0.0,c  )  with  strength  A  in  both  x  and  z  directions 
with  moments  about  y  axis  in  semi  infinite  solid,  the  Galerkin  vector  is  (  Min  dim  and 
Cheng,  1950A) 


-  ♦  x3-c  x3-3c  2c  (x3  +  c) 
F  =  i  [— - „ - S - 


-  4(l-v)(l-2v)  log(R2  +  x3  +  c)] 


-  X[_X[  8v(l-v)xt  4(l-2v){(l-v)x3  -  vc}xt  ^  2cx3xt 

R1  R2  R2  R2^R2  +  x3+c)  d  3 

K2 

2c(x3  +  c)x !  4(l-v)xli 

+  -3  r7~1-  (20) 

K2 

2  2222  222 
Where  Ri  —  (  xt  —  c  )  +x2+x3  and  R2  =  (  Xj  +  c  )  +  x2  +  x3 .  The  stress  field  of  these 

double  forces  can  be  obtained  accordingly  and  by  substituting  Eq.  (18),  the  dispalcement  in 

T  T 

the  matrix  due  to  the  inclusion  with  pure  shear  eigenstrain  e3 ;  =  e  ( 3  .and  all  others  equal  to 
0  is  found  to  be 


e  3  i  in  n 

ui  = - [  V  u3-2(1-v)$  3  ~ZzV  1133“  (3— 4v)y  n3 

4rc(l-v) 

,  n  n  n 
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4n(l-v) 
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u3  = 


T 

e13 
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4jc(  l-V) 


[  y  133-  2(l-v)^tl-  2zy  1333+  (3-4v)y  133 


+  2z  <|>>i33+4vz(j> >13—  2(1— v)<|» tl  ]. 


The  corresponding  stress  field  is 

T 
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and  the  dilatational  field  is 
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(23) 


(l+v)ue13  1  „  n  11  n 

o=  -  — [ - 4*13“  3<(>  13-2z^>  133+2\|r  l333J. 

Jl(l-V) 

T  T 

The  solution  for  inclusion  with  e23  =  e32,  all  others  =  0  is  found  by  cyclic  permutation  of 
(  1,2,3)-  Solution  Eqs.(ll)  and  (22)  can  be  shown  to  satisfy  the  boundary  conditions,  Eq. 
(1) ,  equilibrium,  Eq.  (2)  and  compatibility,  Eq.  (3). 

3.  DISCUSSION 

In  the  present  analysis,  the  elastic  field  caused  by  an  ellipsoidal  inclusion  was 
investigated  for  a  semi-infinite  solid.  An  infinite,  isotropic  elastic  space  was  considered 
first  with  two  ellipsoidal  domains  of  the  same  shape  with  centers  located  at  (  0,0,c  )  and  ( 
0,0,-c  ).  The  two  domains  are  arranged  to  be  mirror  images  of  each  other  with  initial 
eigenstrain  equal  in  magnitude  but  of  opposite  sign.  The  stress  field  for  given  initial 
eigenstrains  in  the  two  domains  can  be  obtained  from  superposition  of  stresses  of  each 
single  domain  in  an  infinite  solid  obtained  by  Eshelby's  method.  A  fictitious  stress  is 
superimposed  such  that  the  the  plane  x3  =  0  is  stress-free.  This  fictitious  stress  is  obtained 

by  the  method  of  Hankie's  transformations  or  by  integration  of  the  elastic  solutions  for 
nuclei  of  strain  in  the  semi-infinite  solid. 

An  important  aspect  that  should  be  pointed  out  here  is  that  the  image  field  ( the  stress 
field  due  to  the  image  inclusion  plus  the  fictitious  stress  )  which  is  superposed  to  the  stress 
field  from  the  single  inclusion  in  the  infinite  solid  in  order  to  satisfy  the  free  surface  free  of 
stress  will  change  the  initial  eigenstrain  in  the  inclusion.  More  specifically,  the  problem 
considered  here  was  that  when  a  initial  eigenstrain  of  the  inclusion  in  an  infinite  solid  is 
given  instead  of  a  given  initial  condition  in  a  half-space  under  consideration.  Therefore  the 
statements  given  in  sections  §  2.1  and  §  4.1  in  Eshelby's  paper  (  1961  )  regarding  the 


sum  of  the  original  field  and  the  image  field  are  chosen  so  that  the  boundary  conditions  on 
the  outer  surface  of  the  matrix  are  satisfied  should  be  carefully  interpreted.  For  example, 
consider  a  cavity  in  an  half-space  under  external  load  as  a  special  case  of  the  ellipsoidal 
inhomogeneity.  Thai  the  image  stresses  not  only  have  to  satisfy  the  boundary  condition  at 
the  free  surface  but  also  vanish  at  the  inclusion-matrix  interface.  Consequently,  when  an 
inclusion  in  the  half-space  is  considered,  the  image  stress  should  satisfy  not  only  the  free 
surface  boundary  conditions,  but  also  the  boundary  conditions  at  the  inclusion-matrix 
interface. 

It  was  also  found  out  by  using  the  relation  I0L=  a/ 2R  when  a  — »0  (  as  given  in 
Eq.(14) ),  we  can  obtain  the  elastic  solutions  of  some  of  the  nuclei  of  strain  in  the  semi¬ 
infinite  solid  by  directely  application  of  Hankel's  transformation  technique  and  the  Kelvin's 
solution  for  a  force  in  an  infinite  body.  Instead  as  shown  by  Mindlin  and  Cheng  (  19S0A  ) 
that  they  are  derived,  by  the  processes  of  superposition,  differentiation  and  integration, 
from  the  solution  for  the  single  force  in  the  interior  of  the  semi-infinite  solid.  The  nuclei  of 
strain  solved  by  present  approach  are  those  axisymmetric  to  the  x3  axis  which  is  normal  to 
the  free  surface,  such  as  single  force  in  x3  direction,  double  force  in  x3  direction,  center  of 
dilatation,  doublet  with  axis  parallel  to  x3  axis. 
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Fig.  1  -  Ellipsoidal  inclusion  Qj  with  principal  half-axes  a,,  a,  and 
in  a  half  space  and  its  image  flj. 
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